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Abstract 

We investigated the kinematic and excitation structure of the NGC 1068 narrow-Hne region (NLR). We 
obtained profiles of several emission lines, [OlIl]A5007, H/?, [Ol]A6300 and [FeVIl]A6087 at high- velocity 
resolution (R ^ 7500 - 11000), and confirmed that they showed different profiles. These profiles are useful 
for understanding the NLR structure, as they cover a wide ionization potential range. By comparing the 
results with a photoionization model, we found that 1) blucshiftcd components at the center are very dense, 
2) those in the northeast region have slightly lower densities than those in the center, and 3) ionization 
parameters of the blucshiftcd components increase with increasing velocity with respect to the systemic 
velocity. We investigated the NLR structure in NGC 1068 based on these results. We show that both the 
observed velocity dependence of the ionization parameter and the gradually increasing velocity field can be 
reproduced by varying the ionizing continuum attenuation, assuming a hollowed biconical geometry and 
varying the column densities of outflowing clouds. 
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1. Introduction 

Narrow-line region (NLR) structure is one of the most 
interesting aspects of Seyfert galaxies and, consequently, 
has been long studied. One tool in such investigations 
is an cmission-linc profile. A blucward asymmetry in 
narrow-line profiles is frequently observed in Seyfert galax- 
ies, and this has led to two traditional interpretations: an 
outflow model including extinction between line-emitting 
clouds, and an inflow model in which extinction originates 
in individual clouds (e.g., Whittle 1985a; De Robertis & 
Shaw 1990; Veilleux 1991a). In some Seyfert galaxies, 
emission lines with higher critical electron density and/or 
higher ionization potential have lower central velocities 
and/or wider line widths. These features have been inter- 
preted as implying that, in the NLR, velocity decreases 
with increasing distance from the nucleus, and that the 
regions closer to the nucleus are denser and in a higher 
excitation state (e.g.. Whittle 1985b; De Robertis & Shaw 
1990; Veilleux 1991b). 

NGC 1068 is one of the most famous Seyfert galax- 
ies and has been studied extensively. The distance to 
NGC 1068 is 14.4 Mpc (1" corresponds to 70 pc.) and 
its systemic velocity is 1148 km s^^ (Bland-Hawthorn et 
al. 1997). Although NGC 1068 has been considered a 
prototypical Seyfert 2 galaxy, its forbidden lines are much 
wider than those of typical Seyfert 2 galaxies, and the 
emission-line profiles show multiple components (Pelat & 
Alloin 1980; AUoin ct al. 1983; Mcaburn & Pedlar 1986). 
Walker (1968) suggested outflow in the NLR of NGC 1068, 
and Cecil, Bland and TuUy (1990) and Arribas, Mcdiavilla 
& Garcfa-Lorenzo (1996) argued that the velocity field of 



the NLR can be explained as a biconical outflow, based 
on tri-dimensional spectroscopic studies. 

The excellent spatial resolution of the Hubble Space 
Telescope (HST) revealed that the striking velocity 
structure of the NLR gradually increased with radius 
(Crenshaw & Kraemer 2000, hereafter CKOO). It must 
be noted that such a velocity field differs from the clas- 
sical perspective, which states that higher velocity com- 
ponents lie closer to the nucleus. The velocity field was 
shown to be reproducible with a hollowed biconical out- 
flow model (Das et al. 2006; Cecil et al. 2002; CKOO), but 
the radial dependence of the velocity in the model was 
included without a physical rationale. Das, Crenshaw & 
Kraemer (2007) and Everett & Murray (2007) attempted 
to explain the radial dependence of the velocity, but no 
successful explanation has been suggested to date. 

The observed emission-line intensity ratios of the 
NGC 1068 NLR can be reproduced by photoionization 
models (Groves et al. 2004; Kraemer & Crenshaw 2000a; 
Kraemer & Crenshaw 2000b, hereafter KCOOb; Alexander 
ct al. 2000; Lutz ct al. 2000). KCOOb decomposed the 
profiles into two kinematic components and investigated 
the physical condition of each component. They reported 
that, while blueshifted clouds suffer from direct ionizing 
continuum, rcdshiftcd clouds are irradiated by an atten- 
uated ionizing continuum. However, the coupling of ex- 
citation structure to kinematics has not been studied in 
detail. 

In this study, wc investigated the kinematic and ex- 
citation structure of the NGC 1068 NLR in detail. We 
performed mcdimii-rcsolution spectroscopy for NGC 1068 
and obtained high velocity resolution [OlIl]A5007, H/3, 
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Fig. 1. FOC image of [O lll]AA4959,5007 (Macchetto et al. 1994; Data were downloaded from the Multimission Archive at the Space 
Telescope Science Institute^.). Two white vertical lines indicate the l."2 slit width aligned for the typical PA=40°. Extraction bins 
(3."4) are also shown. Letters from A to K indicate the cloud locations reported by Evans et al. (1991) and Cecil et al. (2002). 



Table 1. Observation log 



Date 


Object 


Wavelength range 


Included lines 


Exposure time 


PA 




30 Nov 2005 


NGC 1068 


4750 - 5190 A 


H/?, [Olll] 


1800 sec 


37 - 


41° 




NGC 1068 


4750 - 5190 A 


H/3, [Olll] 


1800 sec 


38 - 


39° 




HD 30739 


4750 - 5190 A 




100 sec 






24 Jan 2006 


NGC 1068 


6020 - 6440 A 


[FeVll], [Ol] 


1800 sec 


36 - 


45° 




NGC 1068 


6020 - 6440 A 


[FeVll], [Ol] 


1800 sec 


32 - 


39° 




HD 74280 


6020 - 6440 A 




150 sec 
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[Ol]A6300 and [FeVIl]A6087 emission-line profiles. As the 
ionization potentials of these emission lines encompass a 
wide range, eV to 100 eV, we can research the physi- 
cal condition of the NLR over a wide range of excitation 
states. 

The outline of this paper is as follows: We describe 
observations and data reduction in section 2, and summa- 
rize features of the obtained profiles and velocity depen- 
dences of the emission-line ratios in section 3. In section 4, 
we compare the observational results to a photoionization 
model. In section 5, we discuss the relationship between 
ionization potentials and the central velocities of emission 
lines, and then investigate the kinematic and excitation 
structure of the NGC 1068 NLR. 

2. Observations and data reduction 

2.1. Observations 

We performed medium-resolution spectroscopy of the 
nuclear region of NGC 1068 with a new optical spectro- 
graph (Ozaki & Tokimasa 2005) mounted at the Nasmyth 
focus (F/12) of the 2-m NAYUTA telescope at Nishi- 
Harima Astronomical Observatory on 30 November 2005 
and 24 January 2006. The observation log is summarized 
in table 1. Weather conditions were clear and stable for 
the first night, and almost clear with occasional clouds 
for the second. The seeing was approximately 2 arcsec for 
each day. On 24 January 2006, however, the spatial reso- 
lution degraded significantly because of seeing variations 
or telescope tracking errors during the exposures. 

Slit width and length were 1"2 and 5', respectively. A 
grating with 1800 grooves mm~^ was used. This combi- 
nation of slit and grating provided spectral resolutions of 
7500 at 5000 A, and 11000 at 6220 A, corresponding to 
velocity resolutions of 40 km s^^ and 27 km s^^, respec- 
tively. 

We used the NGC 1068 nucleus as a marker to guide 
the telescope. Unfortunately, the image rotator did not 
function well; hence, the slit swept out some area during 
each exposure due to image rotation. The position angle 
(PA) range for each exposure is listed in table 1. After 
each exposure, we tried to bring the PA back to the ini- 
tial position. We were successful on 30 November 2005 
but not on 24 January 2006. The maximum difference in 
PA was 13°. The NE and SW extraction windows were lo- 
cated at 3"4 from the center; 13° corresponds to 0"77 for 
the NE and SW extraction windows. This difference does 
not affect the results significantly, as the seeing size was 
approximately 2". The typical PA was 40°, which is ap- 
proximately aligned to a radio jet axis (PA=32°; Wilson & 
Ulvestad 1983). Figure 1 shows the average slit position. 

We also observed photometric standard stars HD 30739 
and HD 74280 for flux calibration on the same nights. 



http:/ /archive. stsci.edu/ 
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Fig. 2. Comparison of continuum radial profiles normalized 
by peak intensity. 

2.2. Data reduction 

We reduced the data with IRAF^. We performed bias 
subtraction by referring to overscan regions, and then 
made flat-field corrections using flat-field frames made 
by combining three flat frames obtained after each set of 
exposures. We performed wavelength calibrations using 
comparison frames from a hollow cathode lamp (FcNcAr). 
The root mean square error for the wavelength calibra- 
tion was less than 0.1 A. We checked the systematic er- 
ror of the wavelength calibration with a night-sky line 
of [Ol]A6300.3, and measured the centers of the line as 
6300.1 A and 6300.2 A for each frame. We concluded 
that wavelength calibration accuracy was better than 0.2 
A. After background subtraction, we ran a flux calibra- 
tion and rejected cosmic ray events. Finally, we applied 
the correction for heliocentric velocity. 

Continuum radial proflles around the center of the 
[OlIl]-|-II/3 frame and the [FeVIl]-f[Ol] frame were dif- 
ferent. Comparing g-, r- and i-band images of NGC 1068 
downloaded from the Isaac Newton Group Archive"^, we 
confirmed the absence of any considerable variation in 
the continuum radial profile with wavelength. Thus, 
this difference is thought to have been caused by seeing 
variations or telescope tracking errors during the expo- 
sure. We matched the continuum radial profile in the 
[OlIl]-|-II/3 frame to that of the [Fe VIl] -I- [O l] frame us- 
ing the Gaussian convolution along the spatial direction. 
The resultant spatial resolution after the convolution was 
3 "4. The results of this procedure arc shown in figure 2. 
Spectra at the continuum peak, NE and SW were ex- 
tracted using a bin length of 3"4 along the slit. 

Flux at the center is thought to contaminate flux at the 
NE3."4 and SW3"4, since the flux at the center is much 
stronger. Thus, we estimated central fltix contamination 



IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation, 
http : / /casu. ast . cam .ac.uk/casuadc/archives/ingarch 



4 



S. Ozaki 



Table 2. Selected stars for the template stellar spectra 



Spectral type 


Name 


[Fe/H] 


F5III 


HD 171802 


0.10 


GOIII 


HD 039833 


0.04 


G5III 


HD 027022 


0.13 


KOIII 


HD 130322 


0.06 


K5III 


HD 164058 


-0.07 



for the NE3."4 and SW3"4 and subtracted it from each 
spectrum, assuming that the central spectrum was a spa- 
tially extended Gaussian function with a 3f'4 FWHM. 

We applied an atmospheric absorption-line correction 
around [Ol]A6300. Figure 3a shows a spectrum of 
HD 74280 normalized by its continuum level. As this 
star's spectral type is B3V, it has few intrinsic absorp- 
tion lines, so the atmospheric absorption lines are clearly 
distinguishable. We applied an atmospheric absorption- 
line correction for the 6273 - 6327 A range by dividing the 
spectrum of NGC 1068 by the normalized standard star 
spectrum. Figure 3b shows the spectrum of NGC 1068 
before and after the correction. We can see the negligible 
effect for the [Ol]A6300 profile. 

2.3. Stellar component subtraction 

Stellar components should be subtracted from the ob- 
ject spectra to obtain accurate emission-line fluxes. The 
stellar component is a composite of various stellar spec- 
tral types. However, spectral synthesis is difficult, as the 
spectra in this study cover only a short wavelength range. 
Hence, we assumed that the stellar component in each 
spectral range could be represented by a single spectral 
type. 

We selected some stellar spectra from F5HI to K5HI 
from the ELODIE stellar library (Prugniel & Soubiran 
2001), since bulge light is dominated by late-type red gi- 
ant stars. As the metallicity of the stellar component of 
NGC 1068 is uncertain, we selected stars with metallici- 
ties close to the solar value. The selected stars are listed 
in table 2. 

We created template spectra as follows: Spectra in the 
library were masked at cosmic-ray events or telluric ab- 
sorption lines. The masked regions were typically a few 
pixels. We corrected the masked regions with a linear in- 
terpolation, and shifted the spectra of these stars by the 
recession velocity of NGC 1068. Then we convolved the 
spectra with velocity dispersions of 50, 100, 150, 200, 250, 
and 300 km s~^. We created many template spectra with 
various combinations of spectral types and velocity disper- 
sions. We then scaled the template spectra by referring to 
line-free regions and subtracted the scaled template spec- 
tra from the object spectra. 

We compared the subtracted spectra by eye, and found 
that, when we adopted the G5III spectrum convolved with 
a 150 km s""'^ velocity dispersion, the stellar-origin absorp- 
tion lines in the object spectra were canceled out most pre- 
cisely in both wavelength regions. Therefore, we adopted 



Table 3. Adopted wavelength. 



Line name 


Wavelength 


H/3 


4861.3 A 


[Oni]A5007 


5006.8 A 


[FeVll]A6087 


6087.0 A 


[Ol]A6300 


6300.3 A 



this case as a best-fit template for the underlying stellar 
spectrum. Using GOIII or KOIII with a dispersion of 100 
km s~^ or 200 km s~^ resulted in no significant difference 
from the best-fit template. 

After the template spectrum subtraction, some 
residuals appeared (figure 4), especially around 
[O l]AA6300, 6364. These residuals can be explained 
to some extent by K5III star contamination. We fit these 
residuals with a fourth-order polynomial for each spec- 
trum and subtracted it. The results of this continuum 
subtraction are shown in figure 4. Deviations around 
4750 A and 5170 A are probably due to faint [Arlv]A4740 
and [FeVIl]A5159 emission lines. 

2.4-. Scaling between different wavelength regions 

We scaled the spectra obtained on the different nights 
to reproduce the reddening-corrected [Ol]/[Olll] ratio re- 
ported by Ho, Filippenko & Sargent (1997). In this proce- 
dure, reddening corrections were applied simultaneously. 
Ho, Filippenko & Sargent (1997) adopted the color ex- 
cess E{B — V)= 0.54 derived from the Balmer decrement. 
Koski (1978) reported E{B -V) ^ 0.52 for NGC 1068. 
KCOOb reported E{B -V)= 0.13 - 0.55 and noted that 
the blueshifted components were more heavily reddened 
than the redshifted ones. Thus, the extinction for the red- 
shifted components might be overestimated in this study. 

3. Results 

3.1. Line profiles 

NGC 1068 has been studied extensively and many kinds 
of observations have been performed. However, a detailed 
comparison of several emission lines at high-velocity res- 
olution, such as presented in this study, has not been 
performed except for [Olll] and H/3 (Alloin et al. 1983). 
Figure 6 shows the emission-line profiles. In figure 7, we 
compare the profiles at each position. The [Fe Vll] profiles 
were smoothed with a 3-pixel running mean in figures 6 
and 7. In figure 7, the [FeVll] profile for the SW 3."4 is 
not shown because of its low signal-to-noise ratio. In the 
following subsubsections, we describe the features of the 
line profiles at each position. 
3.1.1. Line profiles at the center 

The [O 111] A5007 profile has four peaks and many shoul- 
ders. These features are consistent with previous high- 
dispersion spectroscopic studies (Meaburn & Pedlar 1986; 
Dietrich & Wagner 1998). The H/3 profile is globally sim- 
ilar to the [O 111] profile. 

An [O l] profile peak was found at the same velocity 
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Fig. 3. Left panel: Normalized HD 74280 spectrum. Right panel: Dashed and solid lines show the spectra at the center around 
[Ol]A6300 before and after the atmospheric absorption line correction, respectively. The vertical axis represents relative strength. 
Flux integration ranges are also shown (subsection 3.2). i 




Fig. 4. Results of the stellar continuum subtraction at the center. Left panels show the object spectra before continuum subtraction 
(dotted line), and the adopted template spectrum (solid line), composed of the stellar spectrum (dashed line) and the polynomial 
(dot-dashed line). Right panels show the spectrum after subtraction. 
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Fig. 5. Same as figure 4, but for the NE3f'4. 



-170 km s"^) as [Olll] and H/3. This means that the 
peak of each line originates in the same cloud and that 
the cloud is radiation-bounded because it has a partially 
ionized zone. 

The [O l] profile displays a narrower shape than the oth- 
ers. In the [O l] profile, we found no well-defined peaks 
around —410 km and —900 km s^^, where [Olll] 
and H/3 exhibit peaks. However we recognized a shoulder 
around —410 km s~^ and a weak tail around —900 km s~^. 
The 4-320 km s^^ peak of the [Ol] profile is much weaker 
than the others. These features suggest that [Ol] flux is 
suppressed in the high- velocity clouds. 

The -1-320 km s~^ peak of the [Ol] profile may be pro- 
duced by [SlIl]A6312 which is separated from [Ol]A6300 
by about -1-560 km s~^. If [Slll] has the same velocity 
as H/3, [Olll] and [Ol], the expected [Slll] peak would 
be located at -1-390 km s~^ with respect to the systemic 
velocity of [Ol]. Its position is marked in figure 6c. A 
displacement between the expected peak location and the 
-1-320 km s~^ peak is clear. No feature was found at the 
expected location. The -1-320 km s~^ peak of [Ol] is con- 
sistent with those of the other lines (figure 7b). Hence, we 
conclude that the contribution from [S III] is negligible. 

The double-peaked [Fe VIl] profile is quite different from 



the other profiles at first glance, as reported by Rodrfguez- 
Ardila et al. (2006). Capetti, Axon & Macchetto (1997), 
in a narrow-band imaging study with HST, showed that 
[NeV] originates from the [O IIl]-bright NLR clouds. The 
ionization potentials of Ne'*"*" and Ne^+ are 97 eV and 
126 eV, respectively, similar to those of Fe^+ (100 eV) 
and Fe^+ (128 eV). Thus, it is plausible that [FeVIl] is 
emitted from the [O III] -bright NLR clouds. 

A detailed comparison of the [FeVIl] profile with the 
others reveals the following three features: First, the 
blueshifted and redshifted [Fe VIl] peaks correspond to the 
[Olll] peaks at —410 km s~^ and 4-320 km s~^, respec- 
tively. Second, no counterpart appears in the [Fe VIl] pro- 
file for the —170 km s^^ peak seen in the other lines. 
Third, in the wavelength region bluer than —410 km s~^, 
[FeVIl] is relatively stronger than the other lines. 
3.1.2. Line profiles at NE 

The [Olll] profile has wide bumps in the bluer part, 
while the [O l] profile declines smoothly. The [O l] profile 
does not exhibit a peak around -1-320 km s~^; the other 
fines do. These features indicate that [Ol] is suppressed 
in high-velocity components like the center. 

A bright narrow component around the systemic veloc- 
ity called a "velocity spike" by Meaburn & Pedlar (1986) 
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Fig. 6. Emission lino profiles. The [Fcvn] profiles are smoothed with a 3-pixel running mean. The horizontal axis is the velocity 
with respect to the systemic velocity, and the vertical axis is relative intensity. Integration ranges are shown in each panel. The 
expected velocity of [Slll]A6312 is shown as a vertical line in the bottom-left panel. 



is predominant except for the [Fe VIl] profile. This compo- 
nent was marginally detected in the [Fe VIl] profile, but its 
relative strength was considerably lower than that of the 
other emission lines. The blueshiftcd peak of the [FcVIl] 
profile is bluer at the NE3"4 than the center, as pointed 
out by Rodriguez- Ardila et al. (2006). 
3.1.3. Line profiles at SW 

The profiles, except [Fe VIl], show similar shapes. [Olll] 
is relatively stronger in the redshifted part than the other 
lines. This result is consistent with KCOOb, who re- 
ported that, in the SW region of the ionization cone, the 
redshifted components show higher excitation than the 
blueshifted components. 

The profiles, except [Olll], have low signal-to-noise ra- 
tios. Therefore, we did not use the SW data in the follow- 
ing sections. 

3.2. Line intensity ratios 

To improve the signal-to-noise ratio, we integrated the 
flux over velocity bins, which were set to include the char- 
acteristic features. Those features and the adopted veloc- 
ity ranges are given in figure 8 and table 4. Hereafter, 



Table 4. Velocity ranges for integration. 



ID 


Feature name 


Velocity range 




1 


—900 km s^^ peak 


-1200 - -800 km s" 


1 


2 


—660 km s^^ shoulder 


-800 - -580 km s" 


1 


3 


—410 km s^^ peak 


-580 - -300 km s" 


1 


4 


— 170 km s^^ peak 


-300 - -50 km s"! 




5 


Systemic 


-50 - +200 km s" 


1 


6 


+320 km s-i peak 


+200 - +400 km s" 


1 


7 


+450 km s~^ shoulder 


+400 - +600 km s^" 


1 


8 


+640 km s^^ shoulder 


+600 - +800 km s" 


1 



the velocity bins are identified with letters indicating the 
location and ID mmrber of the velocity bin, e.g., "CI", 
"NE4" . 

Figure 9 shows the variation in emission-line intensity 
ratios with velocity. Error bars indicate the one-sigma 
level for read-out noise and photon noise. 

Uncertainties in the spectral type and velocity disper- 
sion of the template spectra used in the stellar component 
subtraction were considered. We confirmed that the un- 



S. Ozaki 



[Vol. , 




2000 



2000 



2000 



Velocity (km s' 



Fig. 7. Comparison of the line profiles at each location. Each 
line profile is normalized by its peak intensity. The [Fevil] 
profiles are smoothed with a 3-pixel running mean. The 
[Fevil] profile at the SW3"4 is not shown because of its low 
S/N ratio. Integration ranges arc shown in each panel. 
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Velocity (km s' ) 



Fig. 8. The [OIII] profile at the center. The features referred 
to in the text, the integration ranges and the ID numbers of 
the velocity bins are shown. 

certainties of almost all data points were smaller than, or 
comparable to, the error bars in figure 9 and do not affect 
the results significantly. Two exceptions are CI and NEl. 
When we adopted KOIII as the template spectrum, the 
[Ol] fluxes of CI and NEl were approximately 1.3 times 
larger than the other cases. Since the template spectra 
were shifted by the recession velocity of NGC 1068, those 
spectra included the atmospheric absorption band in the 
CI and NEl velocity range. A stellar-origin absorption 
line is seen between atmospheric absorption lines in the 
KOIII spectrum, which overlaps the atmospheric absorp- 
tion lines in the other templates because of the difference 
in stellar radial velocities. Hence, we may have underesti- 
mated the [Ol] flux of CI and NEl in the adopted G5III 
case. However, even if the [Ol] flux of CI and NEl is in- 
creased to 1.3 times larger than the measured value, our 
discussion is not significantly affected. 

To evaluate contamination of the red wing of 
[OlIl]A4959 into the blue wing of [OlIl]A5007, we shifted 
the [OlIl]A5007 profile by a separation between two lines 
and scaled it down to match the [OlIl]A4959 profile. As a 
result, we confirmed that such contamination is negligible. 

As shown by KCOOb and Cecil et al. (2002), the 
[OlIl]/H/3 ratio is larger in the blueshifted bins than the 
redshifted bins, and this feature is more remarkable at 
the NE3."4 than the center. This was also seen in the 
[FeVIl]/H/3 ratios. The [Ol]/H/3 ratio was largest at C4 
and NE4, and its value monotonically decreased with sep- 
aration from C4 and NE4. 

Figure 10 shows line-intensity-ratio diagrams. In fig- 
ure 10a, the blueshifted components of the center and 
NE3."4 shift to the upper right as velocity increases with 
respect to the systemic velocity. The redshifted and sys- 
temic components of the center and NE3"4 occupy the 
bottom-left region. In figure 10b, the blueshifted compo- 
nents of both the center and NE3"4 shift to the lower right 



NLR Structure of NGC 1068 



9 




1000 



1000 



1000 



Velocity (km s' ) 



Fig. 9. Variation in emission-line intensity ratios with veloc- 
ity. Filled circles connected by a solid line indicate the velocity 
bins at the center, and open circles connected by a dashed line 
indicate the velocity bins at the NE3"4. SW components are 
not shown because of their low S/N ratio. Each data point is 
labeled with letters indicating the position and velocity bin 
listed in table 4. 



with increasing velocity. The data points of the redshifted 
components shift lower with increasing velocity with re- 
spect to the systemic velocity. 

4. Photoionization model 

In this section, we compare the observed results with a 
photoionization model in order to understand the physi- 
cal conditions of each velocity component. Cloudy version 
06.02d (Ferland et al. 1998) was used for the model cal- 
culations. 

Assumptions and input parameters 

We adopted a power law for the ionizing continuum 
shape. Pier et al. (1994) reported that the spectral in- 
dex (a; (X v") of the nuclear ionizing continuum of 
NGC 1068 is —1.7 from interpolation between the UV 
and the soft X-ray data, while KCOOb succeeded in repro- 
ducing the observed emission-line intensity ratios with a 
shallower slope, a = — 1.4. Given such uncertainties, we 
varied a in our photoionization model from —1.7 to —1.4 
with an interval of 0.1. In addition, we varied hydrogen 
density (tih) from 10^ cm~^ to 10^ cm~'^ and the ioniza- 
tion parameter ([/) from 10^"'^ to 10"'^. U is given by: 



U- 



(1) 



47r c rifj ' 

where Q is the number of ionizing photons emitted per 
unit time, r is the distance from the nucleus, and c is 
the speed of light. We assumed a dust-free solar chemi- 
cal composition, as we could not constrain many physical 
parameters with only four emission lines. We assumed a 
constant-density cloud with plane-parallel geometry. The 
termination criterion for the calculation was an electron 
temperature of 1,000 K for enough neutral region to be 
included; i.e., the radiation-bounded case. 

4-2. Results of the model calculation 

In figure 11, the results of the model calculation for 
a = —1.5 are plotted together with the observed data. 
The [Olll]/H/3 ratio reversed at a density of about 10^ 
cm~^ because of [OlIl]A5007 coUisional de-excitation. 
The [Fe Vll]A6087/II/3 ratio reached a maximum at U ^ 
10~^'^ because the He+ absorption suppressed the FeVll- 
emitting region when U > 10^^'^. Owing to this feature, 
this ratio is a good indicator for the ionization parame- 
ter around U = 10~^'^. The [O l]/H/3 ratio increased with 
spectral index, since a partially ionized region becomes 
relatively wide compared with a fully ionized region with 
increasing spectral index. The [Ol]/II/3 ratio decreased 
as the ionization parameter increased. This is explained 
as follows: when the spectral index and density are fixed, 
the size of the partially ionized region does not depend on 
the number of ionizing photons because it approximates 
a mean free path for the ionizing photons, while the size 
of the fully ionized region increases with the number of 
ionizing photons. 

[CaV]A6087 may overlap with [Fe Vll]A6087. To esti- 
mate this effect, we reproduced the plot from figures 11a 
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Fig. 10. Line ratio diagrams, (a) [O III] A5007/H/3 vs. [Fe vil] A6087/H/3. (b) [Olll]A5007/H/3 vs. [O l] A6300/H/3. Filled circles 
connected by a solid line indicate the velocity bins at the center, and open circles connected by a dashed line indicate the velocity 
bins at the NE3f'4. Each data point is labeled with letters indicating the position and velocity bin listed in table 4. 



and lib, except that the vertical axis now represented 
[Fe Vll] + [CaV], and confirmed that [CaV]A6087 had a 
negligible effect on the results. 

4-.3. Comparison with observations 

We adopted a = —1.5 as the best-fit case. When 
a = —1.7 or —1.6, the [Ol]/H/3 ratios predicted by the 
model were smaller than observed. When a = —1.5 
rather than —1.4, the model simultaneously reproduced 
more data points in the [FeVll]/H^ vs. [Olll]/H^ and 
[Ol]/H/3 vs. [Olll]/H/3 diagrams. Parameters derived 
from model/observation comparisons are listed in table 5. 
The comparisons are described in depth in the following 
subsubsections. 

4-. 3.1. The blueshifted components 

C3 and C4 have unique U and nn parameter sets. 

C2 has two possibilities, [U, nn) = (10^^ '^, 10^-^) and 
(10~°-^, 10^-^). When we adopted the higher excita- 
tion case, the estimated luminosity based on those pa- 
rameters was highly inconsistent with the previous result 
(see subsection 4.4). Hence, the lower excitation case, 
([/, tt-h) = (10~^-^, 10^-^), is more plausible. 

Although CI has a unique parameter set, (U , 
nH)=(10~°'^, 10^'^), the estimated luminosity was incon- 
sistent with the previous result, similar to the higher ex- 
citation case for C2. In the [Ol] profile at the center, flux 
contamination from the neighboring bin to CI could be 
relatively large because the feature at CI is very weak. 
If CI is matter-bounded and does not include any [Ol] 
flux, alternative parameter sets for this component can 
be derived only from the [Olll]/H/3 vs. [FcVll]/II/3 di- 
agram: {U, nn) (10"^-^, 10^-^), (10"i-6, lO^'^) and 
(10-11, lO^ *). Given ([/, nn) = (lO^i ^ lO^-^), the esti- 
mated luminosity was roughly consistent with the result 
of Pier et al. (1994) (see subsection 4.4), and we could 
easily explain the difference in central velocity among the 
neutral elements (see subsection 5.1). Thus, wc concluded 
that {U, nn) = (lO"! '^, lO'^-^) was most plausible for CI. 

The blueshifted components at the center are very 



dense (lO^ ^ - 10^-^ cm^^). Walsh (1983) presented the 
[Olll]A4959 and the [Olll]A5007 profiles of the nucleus of 
NGC 1068 at high spectral resolution, and derived densi- 
ties of ^ 10'"' cm~"^ for the blueshifted components, assum- 
ing an electron temperature of 15,000 K. His estimations 
were slightly lower than our results. Since the [O 111] ratio 
depends on both density and temperature, results arc af- 
fected by the assumed temperature. Our results were also 
affected by uncertainty in abundance and dust content. 
In spite of these uncertainties, the difference between the 
derived densities is within a factor of about two. This 
supports our density estimates. 

The ionization parameters of the blueshifted compo- 
nents increase with increasing velocity with respect to the 
systemic velocity. We were able to detect small variations 
due to the strong sensitivity of the [Fe Vll]/H/? ratio to the 
ionization parameter. 

NEl, NE2, and NE3 have slightly lower densities (-- 10^ 
cm"'^) than CI, C2, and C3, although the uncertainty is 
large. These results are consistent with Axon ct al. (1998) 
who reported that the density of cloud G is larger than 
10^-^ cm"'^ based on the [Ariv] ratio. NE4 does not have 
acceptable parameters. 

4.3.2. The redshifted and systemic components 

For the redshifted and systemic components both at 
the center and NE3."4, U - 10""-^ and nn ^ lO^'^ cm'^ 
were obtained from figure 11. However, these results are 
significantly inconsistent with past studies. Walsh (1983) 
derived that component 6 in his paper, which corresponds 
to our C6, has a density of 4.6 x 10^ cm^"^. Pecontal et al. 
(1997) reported that the narrow systemic component of 
the NE region corresponding to our NE5 had a density of 
approximately 400 cm~'^ from the [S ll]AA6717,6731 ratio. 

The reddening correction may affect the results; the red- 
dening for the redshifted components in this study might 
be overcorrected, as we mentioned in subsection 2.4. This 
overcorrection decreases the [FeVll]/H/3 and [Ol]/H/3 ra- 
tios. When we adopted a lower reddening value, the data 
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Table 5. Derived parameters and estimated lower limits of luminosity. 



Bin 


U 


(cm ■^) 


Cloud 
name 


Projected 
distance (pc) 


Luminosity 
(10^^ erg s^i) 


CI 


10-o.y 


lOb.b 


F sub 


56 


28 




10-1-7 t* 


105.3 t* 






2.8 




10-1-6 t 


102.7 t 






0.0087 




10-1-1 t 


102.8 t 






0.035 


C2 


10-1-7 * 


105.3 * 


F 


56 


2.8 




10-0.9 


105.5 






28 


C3 


10-1-8 


105.3 


D 


61 


2.6 


C4 


10-1-9 


105.5 


E 


48 


2.0 


NEl 


10-1-6 


105.0 


H 


126 


8.8 


NE2 


10-1-8 


105.0 


G 


126 


5.6 


NE3 


10-1-9 


105.0 









* Most plausible case (see text). 

t The matter-bounded case (see text). 
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points shifted upward in each diagram of figure 11, which 
meant that the estimated densities did not change signif- 
icantly. Hence, the overcorrection for reddening does not 
significantly affect the density estimation. 

Assuming that high-density clouds with high excitation 
and low-density clouds with low excitation overlap, we can 
reproduce the observed emission-line intensity ratios of 
the redshifted and systemic components (see Appendix). 
It should be noted that this assumption is consistent with 
Walsh (1983) and Pecontal et al. (1997). 

4-4- Luminosity of the central engine 

We derived the luminosity of the ionizing continuum 
(L) as follows: Given a simple power-law continuum with 
a spectral index a, the luminosity per unit frequency in- 
terval can be written = C v°' , where C is a constant. 
Q and L can be expressed by the following forms: 



C ^ 
h V ha 



and 

L = 



Lu dv = 



C 



n + l 



1 



(2) 



(3) 



where is the frequency of the hydrogen ionization limit 
and h is Planck's constant. From the above three equa- 
tions, we find 



L = Att c h ?iH U 



1 



(4) 



Once the distance between the nucleus and the line- 
emitting clouds is determined, we can derive L with U and 
riH estimated in subsection 4.3. Cecil et al. (2002) per- 
formed medium-resolution slit-scan observations for the 
NGC 1068 NLR with HST, and presented the [O III] pro- 
files of each NLR cloud seen in figure 1. On the basis 
of their peak velocities and intensities, we investigated 
the origins of each feature in our profiles. Although 
many clouds appeared in the central extraction window, 
it can be said that CI, C2, C3, and C4 are largely con- 
tributed from clouds F's sub-peak, F, D, and E, respec- 
tively. Clouds G and H and V-shaped filament appear in 
the NE extracted window. Although the profile of the V- 
shaped filament was not presented by Cecil et al. (2002), 
Pecontal et al. (1997) showed that the V-shaped filament 
is redshifted. NEl and NE2 correspond to clouds H and G, 
respectively. We cannot find the corresponding clouds for 
NE3 and NE4. Projected distances between the nucleus 
and each cloud are summarized in table 5. Derived lower 
limits of the luminosity present a wide range as shown in 
table 5. 

When we adopted {U , nH)=(10-i ^ lO^-^) for CI and 
C2, the range of luminosity lower limits moderated. Pier 
et al. (1994) inferred the intrinsic nuclear spectrum of 
NGC 1068 with UV-optical spcctropolarimetric studies 
and X-ray studies, and estimated the ionizing continuum 
luminosity. When their result is recalculated using a dis- 
tance of 14.4 Mpc (they adopted 22 Mpc), the resultant 
luminosity is 



c 



This study 
Marconi et al.(1996) 
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Fig. 12. Dependence of the heliocentric velocity of the emis- 
sion lines on the ionization potential. The filled circles and 
crosses represent the measurements of this study and Marconi 
et al. (1996), respectively, and the dashed line marks the sys- 
temic velocity of NGC 1068. 

Table 6. Neutral line critical densities 




Line name 


Critical density (cm '^) 


[Ol]A6300 


1.5 X 10^ 


[Nl]A5200 


7.0 X 10^ 


[Nl]A5198 


2.2 X 10^ 


[Cl]A9850 


1.4 X 10"^ 



L = 0.36 X 10* 



Jrcfl 
0.01 



D 



14.4 Mpc 



erg s 



(5) 



where the reficction fraction, /lofi, is the ratio of scat- 
tered to intrinsic nuclear light, and D is the distance to 
NGC 1068. Pier et al. (1994) hsted the reflection frac- 
tions derived by several authors; the vahies range from 
0.001 to 0.05. Even if we adopt the minimum reflection 
fraction of 0.001, the luminosities derived from the ([/, 
"■h)=(10~'''^, 10^'^) case for CI and C2 are one order of 
magnitude larger than the above estimation. The lumi- 
nosities derived from the [U , nH) = (10~^'^, 10^ '^) case for 
CI and C2 are consistent with Pier et al. (1994), assuming 
the minimum reflection fraction. 

5. Discussion 

5.1. Dependence of the central velocity on the ionization 
potential 

We verified the dependence of the central velocity of 
the emission lines on the ionization potential shown by 
Marconi et al. (1996). They showed that the higher ex- 
citation lines of the NGC 1068 micleus are more strongly 
blucshifted with respect to the systemic velocity, and that 
the low-excitation lines are located close to the systemic 
velocity. This correlation was also reported by Kraemer 
& Crenshaw (2000a) and Lutz et al. (2000). Figure 12 
shows the correlation between the Gaussian centroid ve- 
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locity and the ionization potential in this study together 
with the results of Marconi et al. (1996). Our data, ex- 
cept [Ol], agree with the results of Marconi et al. (1996) 
within the error limits. 

However, we found a discrepancy in the neutral el- 
ements. Marconi et al. (1996) used [Nl]A5200 and 
[Cl]A9850 for the neutral elements, whereas the neutral 
element in our data was [Ol]. We calculated the criti- 
cal electron densities for these neutral lines with the ionic 
task in IRAF, assuming Te = lO'' K. The calculated densi- 
ties are tabulated in table 6. The critical electron density 
of [O l] A6300 is about two orders of magnitude larger than 
the others. Therefore, coUisional de-excitation may sup- 
press the [Nl] and [Cl] flux, leading to an [Ol] velocity 
that is bluer by comparison with [N l] and [C l] . This inter- 
pretation is consistent with the result that, when we adopt 
(C/, nH)=(10-^ ^ 105-3) for Cl, the blueshifted compo- 
nents at the center have densities higher than 10^ cm~'^. 

5.2. NLR kinematic and excitation structure 

In the previous section, we found that U of the 
blueshifted components increases with increasing veloc- 
ity. Given radiative acceleration, we expected to find that 
U is velocity-dependent. In the radiation-bounded case 
adopted in our photoionization model, the radiative ac- 
celeration can be estimated by 

(6) 



A TT r^ c N nip 

where N is the column density of the cloud and rup is the 
proton mass. From equations (4) and (6), we find 

h vq a U 



a - 



1 I 



(7) 



where / is the cloud depth. This means that, when an 
optically thick cloud is accelerated by radiation from the 
nucleus, its acceleration is proportional to the ionization 
parameter and inversely proportional to cloud depth. We 
assumed the radiation pressure was the driving force of 
the outflow in the following sections, and investigated 
the kinematic and excitation structure of the NLR of 
NGC 1068. 

5.2.1. Ionizing photon attenuation 

The velocity dependence of U is expected when absorb- 
ing matter with varying column densities exists between 
the nucleus and the NLR clouds. In this situation, the 
clouds irradiated by the more attenuated ionizing contin- 
uum are expected to have lower U and lower velocity. 

NGC 1068 probably has ionizing continuum absorbers 
in the vicinity of its nucleus. Alexander et al. (2000) and 
Spinoglio et al. (2005) argued that the observed infrared 
emission-line ratios were reproduced more accurately by 
an attenuated ionizing continuum than by other contin- 
uum shapes. From optical and UV spectroscopic studies, 
KCOOb reported that the redshifted components are ion- 
ized by the attenuated continuum, though the blueshifted 
components are irradiated directly. Their inferred ab- 
sorbers for the redshifted components of NGC 1068 at- 
tenuated the ionizing continuum by about one order of 
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Fig. 13. The estimated velocities, a) Luminosity is varied 
with a fixed launching radius (3 pc) and column density (lO^-^ 
cm"-^). Deprojected radial dependences of the velocity re- 
ported by Das ct al. (2006), Cecil et al. (2002) and CKOO are 
also shown, b) Launching radius is varied with a fixed lumi- 
nosity (3 X 10*^ erg s~^) and column density (10^"^ cm"-^). 
c) Column density is varied with a fixed luminosity (3 x 10^^ 
erg s~^) and launching radius (8 pc). 
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— Blueshifted high excitation clouds 



Blueshifted moderate excitation clouds 




Galactic plane 



Absorbing clouds 



Redshifted clouds 



Fig. 14. Inferred NLR structure of NGC 1068 from the variable attenuation model together with the hollowed biconical outflow 
model. Black and gray arrows show velocities and ionizing radiation whose width means flux, respectively. 



magnitude. Hence, variation in absorbing column den- 
sity seems to be able to account for the difference in U 
observed in NGC 1068. 

To allow for cloud outflow, radiative acceleration must 
overcome gravity. We used the same enclosed central re- 
gion mass distribution for NGC 1068 as Das, Crenshaw & 
Kraemer (2007); it is given by: 



M{r) = 1.5 X 10^ + 7.1 X 10^ r 



-3.2 X 10' 



-2400 



1.5 



(8) 



where r is distance from the nucleus in pc. The first term 
in equation (8) represents the contribution from the su- 
permassive black hole, the second term from the bulge, 
and the last term from the central star cluster. 

When the cloud is assumed to be launched with a ve- 
locity of km s~^ at a distance of tq, we can find the 
velocity at a distance r from equations (6) and (8): 



vir) = ^h{r) + h{r) + h{r) 
where 



^^(^) H 10.3 ^- 1.3) X 10^(1-1 



h{r) 



-6.1 X 10^ log — 



h{r) = -2.2 X 10^ 



2400 



^■0 



ro -f- 2400 



(9) 

(10) 
(11) 
(12) 



where ro is in units of pc, L45 is the ionizing continuum 
luminosity in units of lO''^ erg s"""^, and N2'i is the column 
density of the cloud in units of 10^^ cm~^. When tih = 10^ 
cm~'^ and U = 10~^'^, the column density of the calculated 
cloud is about 10^^ cm~^. However, this value is affected 
by the termination criterion of the calculation. If the cloud 
has more neutral gas, the velocity decreases. It should be 
noted that, if the column density is more than 10^^ cm~^. 



the [O l]-emitting zone is sufficiently included for nn = 10^ 
cm~^ and U ~ 10"^'^. 

We show the estimated velocities given various lumi- 
nosities estimated in subsection 4.4 in figure 13a. The 
estimated velocity range is '^300 km s^^ for luminosities 
derived for the blueshifted components at the center, and 
is smaller than the observed (^ 700 km s~^). This discrep- 
ancy indicates that the observed velocity range cannot be 
reproduced by this simple model. 
5.2.2. Projection effect 

Next, we considered projection effects as well as various 
attenuations of the ionizing continuum. We assumed the 
three-dimensional structure of the NLR in order to resolve 
the projection effect. A hollowed biconical geometry was 
demonstrated for the NGC 1068 NLR (Das et al. 2006; 
Cecil et al. 2002; CKOO). According to that geometry, 
the cone axis is almost perpendicular to the line of sight, 
and the farthest edge of the cone is approximately aligned 
with the galactic disk (see figure 14). 

The hollowed biconical geometry explains the observed 
velocity difference. When the clouds move with nearly 
the same velocity, cloud E, located close to the limb of 
the ionized cone, has a lower velocity than clouds D and 
F, located close to the center of the ionized cone. In the 
hollowed biconical geometry, the blueshifted clouds with 
lower velocity are located nearer to the galactic disk than 
those with higher velocity. If the absorbers lie in the galac- 
tic disk, the absorbing column density decreases outward 
from the galactic disk, resulting in the observed velocity 
dependence of U. 

Figure 13a also shows deprojected radial dependence of 
velocities reported by Das et al. (2006), Cecil et al. (2002) 
and CKOO. We see that, even if we vary the luminosity, 
launching radius, or column density parameters, the grad- 
ual increase in velocity with radius cannot be reproduced 
(figure 13), which is also reported by Das, Crenshaw & 



NLR Structure of NGC 1068 



15 
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Fig. 15. Inferred NLR strueture of NGC 1068 based on the variable column density model. Black arrows show velocities. 



Kraemer (2007). We have to include another perspective 
to explain the velocity field. 
5.2.3. New perspective 

We propose a new perspective to explain the radial de- 
pendence of the velocity. Up to now, we have considered 
the observed velocity field as a time-sequence of cloud 
motion. When clouds with various column densities are 
launched simultaneously at the same radius, clouds with 
lower column densities attain higher velocities and arc lo- 
cated farther from the center for a fixed time (figure 15). 
This situation might give rise to the observed velocity 
field. 

We estimated the velocity field inferred through such 
a situation, and showed the results in figure 16. In this 
model, we adopted the typical luminosity (3 x 10^^ erg 
s^^) estimated in section 4.4. Figures 16a and b show 
the time-dependence of distance and velocity, respectively. 
Figure 16c shows the distances and velocities of clouds 
with various column densities 8 x 10^ years after launch- 
ing. It is impressive that the estimations match the pre- 
viously reported velocity field well. 

There are uncertainties of the luminosity estimations 
in section 4.4 because of our data qualities and uncer- 
tainties of our photoionization model. Those luminosity 
uncertainties affect the results in figure 16. For different 
luminosities, however, we can obtain similar results with 
different launching radii or different time intervals from 
launching. As an example, we also show the results in 
the half luminosity case in figure 16c. It is noted that, in 
extremely low luminosity case, the clouds cannot outflow. 

The above situation could occur when a molecular cloud 
with many cloud cores of various column densities en- 
counters the ionization cone. Ikeda, Sunada & Kitamura 
(2007) observed the Orion-A molecular cloud with the 
H^^CO+( J=l-0) molecular line, which has a high critical 
density (8 x 10^ cm~'^). They identified 236 dense cores, 
and presented a mean core radius of 0.14±0.03 pc and a 
mean density of (1.6±1.2)xl0'' cm~'^. From these values, 
we estimated the core column density as 7 x 10^^ cm~^. 
This estimated column density is roughly consistent with 



the adopted values in figure 16. 

The maximum core diameter in Ikeda, Sunada & 
Kitamura (2007) is 0.46 pc, but the cloud sizes are mea- 
sured as about 10 pc in figure 1. This might suggest that 
the NLR clouds are gatherings of small cloudlets. This 
idea is supported by the filling factor of the NLR clouds. 
Cecil et al. (2002) hsted the areas, the [OlIl]A5007 lu- 
minosities and the [O IIl]A5007/II/3 ratios for each NLR 
cloud. We estimated the filling factor to be very small 
(~ 10^^) from those values and the densities derived in 
section 4. In this estimation, we adopted 1.24 x 10~^^ 
erg cm'^ as '^t^ j h p / ngrip where jnp is the II/3 emission 
coefficient, and Uf. and are electron and proton den- 
sities, respectively (Osterbrock 1989). This is the value 
when Case B and an electron temperature of 10,000 K 
are assumed. We assumed tt-o = rip = nn- 

The timcscale during which an outflowing cloud with 
500 km s^^ passes through the NLR (~100 pc) is about 
10'"' years. This is less than the inferred AGN lifetime 
of 10^ - 10* years (Haiman & Hui 2001 and references 
therein). Hence, if the above situation is accurate, we are 
observing a transient phenomenon. The narrow lines of 
NGC 1068 are abnormally broader than the other Seyfert 
galaxies. This abnormality can be explained by the tran- 
sient phenomenon of dense cloud outflow. 

In conclusion, we can reproduce both the velocity de- 
pendence of U and the velocity field of the NLR of 
NGC 1068 through varying ionizing continuum attenu- 
ations, a hollowed biconical geometry, and this new per- 
spective. 

6. Summary 

We performed medium-resolution spectroscopic obser- 
vations of NGC 1068 and obtained emission-line profiles of 
[O III]A5007, H/3, [Fe VII]A6087 and [O l]A6300. At the cen- 
ter, [Qui], H/3 and [Ol] exhibited a peak at —170 km s~^ 
with respect to the systemic velocity. The double-peaked 
profile of [Fe VIl] was confirmed as reported by Rodriguez- 
Ardila et al. (2006). Its two peaks at the center were lo- 
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Fig. 16. Time dependence of distance (a) and velocity (b) of 
the clouds with various column densities, when L45 = 3.0 and 
ro = 8 pc. In the bottom panel (c), solid circles show esti- 
mated distances and velocities of clouds with various column 
densities after 8 X 10* years since launching. Labels with the 
circles show the column density in units of lO'^^ cm~^. Open 
triangles show the results in the half luminosity case in which 
ro = 8 pc and time interval since launching is 1 X 10^ years. 
The radial dependence of the velocity reported by Das et al. 
(2006), Cecil et al. (2002) and CKOO is also shown. 



cated at the sub-peaks of [O III] and H/3. We marginally 
detected the narrow spike component in the [FeVIl] pro- 
file at the NE3."4. In the [Ol] profile at the center, the 
features of the high- velocity components were very weak. 

Comparing observations with a photoionization model, 
we investigated the physical conditions of the line-emitting 
regions. The ionization parameters of the blueshifted com- 
ponents increased with increasing velocity with respect 
to the systemic velocity. The bluest components at the 
center might be matter-bounded. The densities of the 
blueshifted components at the center were ~ lO^''^ - 10^'^ 
cm~'^ and those at the NE3"4 were slightly lower than 
at the center. The systemic and redshifted components 
might be constructed from two or more regions with dif- 
ferent excitation states. 

When we assumed a hollowed biconical geometry in 
which absorbing matter with varying column densities ap- 
pears between the nucleus and the NLR clouds, we suc- 
ceeded in reproducing the velocity dependence of the ion- 
ization parameter. However, this model did not explain 
the gradual increase in the velocity field with radius. We 
succeeded in reproducing the velocity field given clouds 
with various column densities launched simultaneously at 
the same radius. Finally, we show that both the observed 
velocity dependence of the ionization parameter and the 
gradually increasing velocity field can be reproduced by 
varying the ionizing continuum attenuation, assuming a 
hollowed biconical geometry and varying the column den- 
sities of outflowing clouds. 
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Appendix. Clouds located along a line-of-sight 

In this appendix, we show the expected location of ob- 
served data in an emission-line intensity ratio diagram 
given clouds with different physical conditions located 
along a linc-of-sight. For simplicity, we consider two 
clouds, clouds 1 and 2. We write fluxes of H/3 and [Olll] 
as F(RP) and F{[0 III]), respectively. The total fluxes are 
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Fig. 17. Line ratio diagrams. Model grid and observed data points are the same as figures 11a and d. Open and filled squares 
indicate, respectively, the locations of each cloud and data point that may be observed. 



given by 

i^tot(H/3) = Fi(H/?)+F2(H/3) 
and 

^^tot ( [O III] ) = Fi ( [O III] ) + F2 ( [O III] ) 



(Al) 



(A2) 



Assuming the ratio of the H/3 flux of cloud 1 to total H/3 
flux is 77, we get 



Fi(H/?) F^im , 

~ "'^'^ 7^ — , - = 1 — ?7 



(A3) 



Ftot(H/?) ' Ftot(H/3) 
From the above equations, we flnd 

^tot([Oiii]) j^i([Oni]) , , F,{[Om]) 

This equation is also applicable to other emission-line in- 
tensity ratios, [FeVIl]/H/3 and [Ol]/H/3. This result indi- 
cates that, when we connect the locations of each cloud 
in the diagram with a line, the observed data are located 
at a point dividing the line into 77: (1 — rj). 

For example, when {U, nH) = (10^^'^, 10^) for cloud 1 
and (lO""^'^, 10'^) for cloud 2, and when those clouds con- 
tribute equally to the total H/3 flux, we expect that the 
observed locations are the black squares in figure 17. 
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